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Pre  face 


The  study  of  radioactive  fallout  from  nuclear  detonations 
is  of  great  interest  in  this  time  of  debate  over  the  proposed 
MX  missile  field.  A  counterforce  attack  against  this  complex 
could  consist  of  thousands  of  large  yield  nuclear  weapons 
detonated  on  or  near  the  ground.  Such  an  attack  could  have 
more  far-rcaching  consequences  than  most  military  planners 
consider.  A  highly  survivable,  mobile  missile  system  would 
have  little  value  to  the  millions  of  people  killed  from  the 
fallout  produced  by  such  a  counterforcc  attack. 

Presented  within  is  a  simple,  efficient  procedure  for 
accurately  determining  this  collateral  damage  of  fallout  for 
any  scenario  involving  many  bursts.  This  method  is  designed 
for  an  operational  planner  to  easily  "scope  the  problem" 
without  utilizing  much  computer  time. 

I  am  grateful  to  Dr.  Charles  J,  Bridgman  for  his  guidance 
in  the  development  of  this  procedure,  and  to  my  lovely  wife, 
Michaela,  for  her  patience  and  support. 

John  P.  Grand ley,  Jr. 


(This  thesis  was  typed  by  Sharon  A.  Gabriel) 
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Abstract 

A  method  is  developed  for  calculating  fallout  deposition 
downwind  from  a  massive  nuclear  attack  on  a  small  target 
area  over  a  short  time  span.  This  is  accomplished  using 
existing  smear  codes  and  replacing  their  existing  horizontal 
activity  distribution  with  an  approximating  function.  This 
function  is  the  difference  between  two  cumulative  normal 
functions  which  are  shown  to  result  from  superposition  of 
individual  bursts.  A  comparison  is  made  between  the  contours 
predicted  by  this  new  code  and  contours  predicted  by  the  old, 
time-consuming,  iterative  procedures.  The  new  code  has  been 
employed  in  several  different  scenarios  involving  the  proposed 
MX  field  to  determine  the  resulting  dose  contours  from  a 
massive  attack  against  that  field. 
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A  MULTIBURST  FALLOUT  MODEL 


FOR  OPERATIONAL  TYPE  STUDIES 

I .  Int  roduct ion 

A  new  Intercontinental  Ballistic  Missile  (ICBM)  system 
will  be  deployed  in  the  Western  United  States  within  the 
next  several  years.  This  new'  system,  designated  the  MX, 
will  greatly  enhance  the  survivability  of  the  land-based 
leg  of  the  Triad.  This  survivability  will  be  accomplished 
through  mobility;  a  single  missile  will  be  shuttled  ran¬ 
domly  among  23  shelters  on  a  racetrack.  The  current  propo¬ 
sal  calls  for  200  missiles  to  be  purchased,  with  200  race¬ 
tracks  to  be  built  in  Utah  . and  Nevada  (Ref  1:3). 

This  method  of  "hiding"  the  MX  missiles  creates  tremen¬ 
dous  problems  in  targeting  for  any  enemy.  One  possible 
targeting  option  would  be  to  attempt  to  destroy  as  many  MX 
shelters  in  as  short  a  time  as  possible.  This  option  would 
have  a  high  probability  of  destroying  a  large  percentage  of 
the  MX  missile  force,  while  interfering  with  the  command, 
control,  and  communication  necessary  to  effect  a  retaliatory 
strike.  Also  of  great  concern  is  the  corresponding  downwind 
fallout  effects  from  such  a  massive  strike.  There  are  several 
single  burst  nuclear  fallout  models  available  which  could 
predict  this  collateral  effect  by  superposition  of  individual 
burst  fallout  patterns. 


Presently  there  are  many  different  individual  burst 
models  being  used  by  many  different  government  agencies.  Mo 
of  these  models  have  their  analysis  based  on  a  study  done  by 
the  Weapons  Systems  Evaluation  Group  (lVSEG-10)  in  1959. 
However,  Bridgman  and  Bigelow  (Ref  2)  have  recently  shown 
that  the  heart  of  the  WSEG  calculations  is  in  error.  They 
proposed  an  alternative  method  of  calculation,  which  was 
used  by  Colarco  (Ref  3)  to  construct  an  improved  single¬ 
burst  fallout  code  for  operational  type  studies.  Colarco's 
code  produces  results  which  are  closer  approximations  than 
WSEG-10  to  the  Defense  Land  Fallout  Prediction  System 
(DELFIC).  DELFIC  is  considered  by  many  to  be  the  benchmark 
of  fallout  computer  codes.. 

One  of  the  unfortunate  aspects  of  these  single-burst 
codes  is  that  they  cannot  easily  compute  the  effects  of 
many  bursts  within  a  relatively  small  target  area  (such 
as  an  MX  field).  Fallout  analysis  of  a  counterforce  attack 
on  a  missile  field  (or  similarly  distributed  target)  is  now 
done  by  superpositioning  hundreds  of  individual  bursts  by 
an  iterative  procedure.  3ecause  of  the  computer  time  needed 
for  these  iterations,  this  is  a  slow  and  expensive  process. 
Also,  because  of  computer  limitations,  this  may  not  be  a 
completely  accurate  procedure. 

In  the  following  chapters,  a  method  will  be  developed 
to  produce  a  fast-running  computer  code  to  accomplish  multi¬ 
burst  calcualtions  without  the  need  for  superpositioning. 


This  code  will  be  designed  for  operational  use;  that  is, 
it  will  be  a  fast,  inexpensive  tool  for  operational  analysts 
to  use  in  predicting  a  reasonably  accurate  fallout  deposition. 
This  code  will  then  be  applied  to  a  counterforce  scenario 
against  the  MX  field  in  Utah/Nevada.  The  results  from  this 
study  will  then  be  compared  to  the  results  obtained  from 
using  WSEG-10  calculations  in  the  same  scenario. 


1 1 .  Calculation  of  the  Multiburst  Distribution 

US EG- 10  predicts  fallout  dose  rates  on  the  ground  with 
the  following  equation  (Ref  4:17): 

00 

61(x,y)  =  k  J  f  (x  ,y ,  t ' )  g(t')  dt"  (1) 

o 

where  k  is  a  source  normalization  constant,  g(t')  is 
the  activity  deposition  rate,  and  f(x,y,t")  is  the 
normalized  horizontal  activity  function.  This  activity 
function  is  a  bivariate  normal  function  with  a  time  varying 

standard  deviation  in  the  cross  wind  direction: 

x-V  t  v 

,  -«-37-5  , 

f(x,y,t)  =  -  e  •  -= -  e 

/2tt  a  >^2tt  a  ft) 

a  y 

(2) 

where  x  is  the  downwind  distance,  y  is  the  crosswind 
distance,  t  is  time  and  a  is  the  standard  deviation. 
is  assumed  to  be  a  constant  wind  which  rotates  uniformly 
with  altitude.  This  results  in  a  constant  value  of  wind 
shear  for  the  determination  of  a  (t) 

It  can  be  seen  that  Eq  (2)  can  be  written  as 

f(x,y,t)  =  f(x,t)  •  f(y,t)  (3) 
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where  t  is  arrival  time  of  the  cloud  and 
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f  (y,t) 
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(5) 


It  is  obvious  from  F.q  (5)  that  f(y,t)  is  a  normal 
function  describing  the  crosswind  spread  of  a  single  nuclear 
cloud.  If  two  or  more  clouds  arc  in  close  proximity  and 
merge  at  some  point  to  become  one  large  cloud,  a  single 
normal  function  will  no  longer  account  for  horizontal 
activity;  there  will  be  some  cumulative  effect  from  each 
contributing  single-burst  cloud.  The  fallout  on  the  ground 
from  this  new,  large  cloud  will  still  be  described  by 
liq  (4),  except  that  f(y,t)  will  no  longer  be  given  by 
Eq  (5). 

f(y,t)  for  Multiple  bursts 

If  an  observer  were  standing  some  distance  downwind 
from  a  target  complex  (for  example,  a  missile  field)  and 
was  able  to  observe  a  large-scale  nuclear  attack  on  that 
complex,  he  would  initially  sec  many  single  nuclear  clouds 
rising  from  the  surface.  If  this  particular  attack  were 
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1  •  >.;k‘  mi-,.',  at  on  Jet  er.a  t  i  or.s  ,  all  go  in.;  off  wi  thin  a  very 

short  time  of  one  another,  and  distributed  uniformly  along 
an  SO  mile  line  perpendicular  to  the  wind  direction,  then 
the  observer  would  see  100  individual  overlapping  clouds. 

If  the  horizontal  activity  distribution  within  each  cloud 
along  this  line  is  given  by  Fq  (S),  then  the  observer  would 
sec  100  normal  functions,  as  in  Figure  1.  Note  that  y  is 
now  the  crosswind  distance  as  measured  from  the  center  of 
the  target  field. 


It  can  be  seen  from  Figure  1  that  any  point  downwind 
from  this  line  of  clouds  will  receive  activity  from  not 
only  the  single  cloud  directly  upwind  from  this  point, 
but  also  from  adjoining  clouds.  This  additional  activity 
becomes  even  more  pronounced  as  tire  clouds  begin  drifting 
with  the  wind  as  the  horizontal  activity  standard  deviation, 
0^  ,  increases  with  time  (as  per  the  WSEG-10  analysis). 

The  total  activity  at  any  point  downwind  can  be  found  by 
finding  tire  contribution  of  each  single  cloud  and  adding 
all  100  contributions  together.  If  the  standard  deviation, 
0^  ,  is  large  with  respect  to  the  intercloud  distance, 

the  addition  operation  can  be  replaced  by  an  integral: 


f(y,t) 
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whore  w/2  is  the  half- field  width  (in  the  example,  40 
miles)  and  N  is  the  total  number  of  bursts  (here,  100). 
This  function  can  be  easily  evaluated  if  it  is  split  into 
two  part  s  : 
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which  is  easily  transformed  into 
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■  y+w/2 

o~Tt7 

I  >  _L 


’ -iZ  ‘ 


/  2tt 


dz 


y  -  w/  2 

VTTT 

f  y  i  -i<7 2 

J  — —  e  *  d  z 

-00  /  2ll 


(8) 


by  letting 


y-yG 


(9) 


and 


dz 


-dy 
1  o 

OTT 


(10) 


8 


Thus,  f(y,t)  Tor  many  bursts  is  the  difference  be l worn 
the  cumulative  normal  functions  for  two  different  argumcnl s . 
When  this  function  is  used  in  f.q  (4),  the  total  dose  rate 
for  any  point  downwind  can  be  found.  A  graph  of  this 
function  is  presented  in  figure  2,  with  the  100  single 
cloud  distributions  added  for  comparison. 

Computer  solutions  for  dose  rates  from  F.q  (4)  would  be 
quick  and  simple  except  for  the  integrals  in  the  multi - 
burst  f(y,t)  .  Therefore,  the  following  approximation 
was  used  for  the  integrals  in  Fq  (8)  (Ref  9:932): 

P(z)  =  1  -  j  (.1968342  +  .11519 4 z2 

+  .  0  0  0  3  4  4  z  3  H  .019S27z**)-t*  (11) 

where  z  is  the  upper  limit  of  each  integral.  The  difference 
between  the  two  integrals  was  then  multiplied  by  the  number 
of  detonations  and  divided  by  the  field  width  to  produce  a 
numerical  answer  for  f(y,t)  .  This  number  was  then  used 
in  Eq  (4)  to  produce  dose  rates. 

As  will  be  seen  in  the  next  chapter,  dose  rates  from 
Eq  (4)  using  this  new  f(y,t)  agree  very  well  with  dose  rates 
derived  by  addition  or  superposition.  However,  several  cau¬ 
tions  need  to  be  mentioned.  This  f  (y)  function  is  meant  to 
be  used  in  situations  where  there  arc  many  bursts  in  a  rela¬ 
tively  small  physical  area.  Also,  the  individual  bursts 
positioned  on  the  field  width  line  cannot  he  more  than  one 
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standard  deviation  apart;  that  is,  the  number  of  bombs 
detonated  in  a  defined  area  (the  bomb  density)  must  be 


large  enough  to 
approximate  the 
must  be  treated 


allow  the  multiburst  f(y)  function  to 
cumulative  effect.  Smaller  bomb  dens  it i 
as  an  aggregation  of  single  bursts. 
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III.  Validation  of  the  Multiburst  Distribution 

Any  model  must,  by  necessity,  be  shown  to  actually 
conform  to  the  reality  it  purports  to  represent.  There  are 
several  ways  in  which  a  model  can  be  validated,  with  most 
validation  methods  comparing  model  results  to  actual  data. 

As  there  is  little  data  available  on  downwind  fallout 
dosages  from  multiple  bursts,  a  different  approach  was 
necessary.  This  involved  a  comparison  between  the  results 
of  the  new  model  and  results  derived  from  an  iterative 
superpositioning  method  of  the  same  scenario. 

Scenario 

A  scenario  involving  a  counterforce  attack  against  the 
missile  field  at  Whiteman  Air  Force  Base,  Missouri  (Ref  7:32) 
was  used  as  a  basis  for  comparison.  This  particular  missile 
field  encloses  150  silos  in  an  area  roughly  90  miles  on  a 
side.  Two  one-megaton  devices  were  then  simultaneously 
detonated  at  each  silo,  for  a  total  of  300  bursts.  These 
300  bursts  were  then  distributed  evenly  on  a  98  mile  line 
perpendicular  to  the  average  west  wind  of  20  miles  per  hour, 
and  a  shear  of  one  per  hour.  Fifty  percent  of  the  yield  of 
each  burst  was  from  fission.  The  coordinates  of  the  contour 
lines  of  the  following  unit  time  reference  dose  rates  were 
then  computed:  1000  roentgen/hour  (r/hr),  500  r/hr,  and 
100  r/hr. 
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Description  of  Superpositioning  Procedure 

A  single  burst  code  similar  to  Colarco’s  model  (Ref  3) 
was  used  for  the  iterative  procedure.  A  grid  was  established 
with  2000  miles  on  the  x  (downwind)  axis  and  182  miles  on 
the  y  (crosswind)  axis.  The  x  axis  was  broken  down  into 
201  increments,  while  the  y  axis  was  divided  into  27. 

For  more  efficient  computer  operation,  the  150  silos  were 
broken  down  into  15  groups  of  10  silos  per  group.  These 
15  groups  were  then  placed  on  grid  line  x  =  6  and  on  every 
grid  line  between  y  =  7  and  y  =  21  ,  resulting  in  15 

silo  groups  evenly  distributed  on  a  north-south  line  98 
miles  long.  To  simulate  20  bursts  at  each  silo  group,  the 
computed  single-burst  dose  rate  at  each  silo  group  was 
multiplied  by  20.  The  superpositioning  method,  then,  was 
accomplished  by  stepping  from  point-to-point  on  the  grid 
and  calculating  and  adding  the  contribution  of  dose  rate 
from  each  silo  group  to  that  point. 

Machine  limitations  dictated  the  coarseness  of  the  grid. 
Consequently,  there  was  a  reduction  in  accuracy.  A  summation 
of  all  grid  points  did  not  total  up  to  the  entire  activity 
produced  by  the  detonations;  approximately  five  percent  of 
the  activity  was  missing.  This  missing  activity  was  not 
deemed  important  to  the  comparisons,  and  is  believed  to  be 
the  result  of  the  coarse  grid  at  the  extreme  downwind  limits 
of  the  fallout  field. 
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Description  of  Multiburst  Procedure 


Results  using  the  new  model  were  obtained  using  the 
multiburst  code  shown  in  Appendix  A.  The  required  inputs 
were  inserted,  and  the  coordinates  for  the  desired  contours 
produced.  No  machine  limitations  were  encountered  using 
this  model. 

Comparison  of  Multiburst  to  Superposition 

The  comparison  was  made  using  the  same  inputs:  300 
one-megaton  bursts,  each  with  50  percent  of  its  yield 
from  fission;  wind  of  270/20  with  a  shear  of  one  hr  ; 
field  width  of  08  miles.  The  contours  obtained  were  for 
doses  of  100,  500,  and  1000  r/hr. 

Two  different  sets  of  contours  were  drawn  for  two 
different  particle  size-activity  distribution  curves.  The 
WSEG  contours  use  a  size-activity  curve  with  parameters  of 
a  =  44  and  £n3  =  .690  .  A  second  set  of  contours  was 

generated  using  size-activity  curve  parameters  of  a  ~  37 
and  £.n 3  =  1.5  28  .  These  parameters  more  closely  resemble 

those  used  by  the  DbLFIC  model  (see  Figure  31.  All  calcu¬ 
lations  using  these  parameters  and  the  Bridgman/Bi gelow 
procedure  for  activity  deposition  rate  (Ref  2:29-30)  will 
be  termed  AFIT  calculations. 

Figures  4a  through  4c  show  the  three  desired  contour 
comparisons  using  the  WSF.G  formulation  for  activity  depo¬ 
sition  rate,  while  Figures  5a  through  5c  show  the  same 
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Figure  4a.  WSFG  Comparison 
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Figure  4c.  WSHG  Comparison 
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comparison  using  the  Al-'IT  activity  deposition  formulation. 

As  can  be  seen,  these  comparisons  of  superpositioning  and 
the  multiburst  code  are  in  good  agreement.  For  the  smaller 
dose  rates  there  is  a  difference  in  the  crosswind  deposition 
of  activity  as  downwind  distance  increases  from  ground  zero. 
This  discrepancy  is  due  to  the  way  the  300  bursts  are  distri¬ 
buted  in  the  superpositioning  procedure;  that  is,  the  15 
groups  of  20  detonations  more  closely  resemble  fifteen 
20-megaton  detonations.  This  aggregate  cloud  will  taper  off 
much  more  than  an  aggregate  cloud  of  300  one-megaton  bursts, 
as  shown  by  Figure  5c.  Consequently,  the  multiburst  contours 
are  considered  to  be  more  accurate. 
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To  demonstrate  the  ease  with  which  this  new  model  can 
be  used,  the  fallout  patterns  from  a  counter  force  attack  on 
an  MX  field  were  generated.  Several  different  scenarios 
we re  devised,  and  the  appropriate  inputs  made  to  the  code 
MULTI,  which  is  listed  in  Appendix  A.  MULTI  incoiporatcs 
the  multiburst  f(y,t)  distribution  developed  in  Chapter 
II.  An  interesting  feature  of  this  cede  is  the  ability  to 
use  either  the  V.’SEG  g(t)  (activity  deposition  rate)  or 
the  AFIT  g(t)  in  dose/dose  rate  computations.  Therefore, 
for  all  scenarios,  a  comparison  study  was  made  of  contours 
generated  using  the  two  different  g(t)'s 

Scenarios 

The  MX  field,  consisting  of  4600  shelters,  is  proposed 
to  be  situated  in  the  states  of  Utah  and  Nevada.  Several 
different  scenarios  were  created  around  these  shelters, 
based  on  several  parameters.  These  parameters  include  the 
number  of  attacking  reentry  vehicles  (RV) ,  the  yield  and 
fission  fraction  of  each  RV ,  and  the  average  continental 
winds.  Two  different  dose  contours  were  generated  for  each 
dose  using  the  two  different  g(t)'s  .  Based  on  all  these 
variables,  a  total  of  64  plots  (or  32  comparisons)  were 
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created . 


To  simulate  a  full,  counterforce  attack,  4600  RV's 
were  targeted  against  the  MX  field.  The  crosswind  width  of 
this  field  was  estimated  at  190  miles.  The  yield  of  each  RV 
was  allowed  to  be  one  of  two  values:  either  one  megaton  or 
500  kilotons.  These  yields  are  commensurate  with  the  war¬ 
head  yields  of  Soviet  ICBM's.  A  50  percent  yield  due  to 
fission  was  assumed  for  each  warhead.  As  another  comparison, 
the  number  of  RV’s  was  decreased  to  2300  with  warhead  yields 
remaining  the  same. 

Several  average  winds  were  obtained  (Ref  8)  for  the 
continental  United  States  as  a  whole  at  an  altitude  of 
40,000  feet.  All  winds  were  from  the  west  (270°)  and  had 
different  velocities  based  on  the  season.  An  average  summer 
wind  had  a  velocity  of  35  miles  per  hour  (or  written  as 
270/35),  while  an  average  winter  wind  was  given  as  270/77. 
Strong  seasonal  winds  were  derived  by  adding  one  standard 
deviation  to  the  average  winds.  This  resulted  in  a  strong 
summer  wind  of  270/70  and  a  strong  winter  wind  of  270/119. 
These  winds  were  then  applied  to  each  scenario. 

Two  different  dose  contours  were  generated  for  each 
scenario,  the  doses  being  1500  and  500  rems.  Since  most 
single  story  residences  above  ground  have  a  protection  factor 
of  three  (Ref  7:33),  the  1500  rem  contour  represents  500  rems 
indoors.  Five  hundred  rems  is  considered  to  be  the  dosage 
necessary  to  produce  50  percent  fatalities  in  30  days 
(Ref  7:32). 
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Fatalities  will  not  be  estimated  in  the  following 
comparisons.  However,  the  contour  lines  do  enclose  the 
areas  of  high  fatalities  from  fallout.  Major  cities  down¬ 
wind  from  the  bursts  are  delineated  on  the  figures  for 
easier  reference.  It  will  be  obvious  that  the  WSEG  formu¬ 
lation  significantly  underestimates  these  areas  of  high 
fatalities . 

Results 

Sixty-four  computer  runs  were  made  interactively 
using  a  CDC  6600  computer.  Average  compilation  time  was 
1.8  seconds  with  average  execution  time  of  1.5  seconds. 

The  64  runs  were  then  combined  into  32  comparative  graphs. 

These  32  figures  were  arranged  in  four  groups  according 
to  wind  velocity.  Figures  6a  through  6h  have  a  common  wind 
of  270/35  (average  summer).  The  second  group,  Figures  7a 
through  7h,  has  a  wind  of  270/77  (average  winter),  while 
Figures  8a  through  8h  have  a  slightly  lower  wind  of  270/70 
(strong  summer).  The  strong  winter  wind  group.  Figures  9a 
through  9h,  is  last  with  a  wind  of  270/119.  All  winds  have 
a  shear  of  1/hr. 

It  is  obvious  from  all  comparisons  that  the  WSEG 
contour  lines  encompass  much  less  area  than  the  AFIT  contour 
lines.  This  is  a  direct  result  of  the  Bridgman/Bigelow 
method  of  computing  g(t)  and  of  the  parameters  used  for 
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the  particle  size-activity  distribution  curve.  Colarco 
showed  the  same  kind  of  result,  only  for  a  single-burst 
model . 


It  is  interesting  to  note  that  the  WSEG  contours  are 
probably  close  approximations  to  the  contours  generated  by 
the  Department  of  Defense  (DoD)  in  their  predictions  of 
fatalities  from  an  MX  field  attack;  the  same  basic  WSEG 
formulation  is  used.  Consequently,  the  estimation  of 
resulting  fatalities  is  not  as  high  as  it  would  be  if  the 
AFIT  contours  were  used,  especially  in  the  low  wind  situa¬ 
tions.  As  the  wind  increases  and  the  WSEG  contours  finally 
stretch  into  the  Atlantic  Ocean,  the  number  of  fatalities 
predicted  by  WSEG  and  AFIT  will  be  about  the  same  as  the 
enclosed  population  is  equal. 
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Figure  6a.  Average  Summer  Kind  Contours 
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Figure  6c.  Average  Summer  Kirn]  Contours 
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Figure  fid.  Average  Summer  Wind  Contours 
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Figure  6e .  Average  Summer  Kind  Contour:'. 
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Figure  6f.  Average  Summer  hind  Contour;- 
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Figure  Oh.  Average  Summer  Wind  Contours 
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Figure  7a.  Average  U'intcr  H'ind  Contours 
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Figure  7b.  Average  V, 'inter  bind  Contour.0 
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Figure  7c.  Average  Winter  Wind  Contours 
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Figure  7d.  Average  Winter  Wind  Contou.s 
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Figure  7c.  Average  Winter  Wind  Contours 
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Figure  7f.  Average  Winter  Kind  Contour:; 
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Figure  7g.  Average  Winter  Wind  Contours 
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Figure  7h.  Average  Winter  Kind  Contours 
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Figure  8a.  Strong  Summer  Wind  Contours 
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Figure  8c. 
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Figure  8d.  Strong  Summer  Wind  Contours 
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Figure  8e.  Strong  V.’ummcr  Kind  Contour 
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Figure  8g.  Strong  Summer  Wind  Contour; 
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Figure  8h.  Strong  Summer  Kind  Contour- 
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rigurc  9a.  Strong  Winter  Wind  Contours 
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Figure  9d .  Strong  Winter  Wind  Contours 
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Figure  9c.  Strong  V.'intcv  Kind  Contours 
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Figure  Of. 


Strong  Winter  Wind  Contours 
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Figure  9g. 


Strong  V.’inter  Wind  Contours 
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Figure  9h.  Strong  Winter  Wind  Contours 
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Recommendat ions 


The  presented  code  will  give  realistic  approximations 
of  the  downwind  deposition  of  fallout  with  great  speed  and 
ease.  This  is  exactly  what  operational  planners  need  to 
make  timely  decisions.  However,  several  enhancements  could 
be  made  to  the  program  to  increase  its  accuracy. 

The  first  enhancement  involves  wind.  The  use  of  a 
single,  constant  wind  from  the  surface  to  40,000  feet  and 
across  the  countryside  is  a  notoriously  bad  approximation. 
Some  way  must  be  devised  to  incorporate  several  different 
winds  into  the  model. 


60 


Hi  h  1  j  o': rapliv 
-  --  - 1  ■■ 


1.  Korb,  Lawrence  J.  ’’The  Case  for  the  MX,"  Air  University 
Review,  31:  2-10  (July -August  1980). 

2.  Bridgman,  C.J.  and  IV.  S.  Bigelow.  "A  New  Fallout 
Prediction  Model  for  Use  in  Operational  Type  Studies," 
Unpublished  draft.  Wright  -  Pa 1 1 ersun  AFB  OH:  Sci'.oo  1 

of  Engineering,  Air  Force  Institute  of  Tech.no  logy . 

3.  Colarco,  Richard  F.  A  Computer  Fallout  Model  for 
Operational  Type  Stud:  es  ~  MS  The  sTs~  W  r  igi  7t  -Patterson 
AFB  01 H  School  oT  1  in  gTTTee  ring ,  Air  };orce  Institute  of 
Technology,  March  1980. 

4.  Pugh,  George  F..  and  Robert  J.  Culiano.  An  Ana 3  yt  ica  1 
Model  of  Close-In  Deposition  of  Fallout  l"7-!  i  s  c  ,i  n 
Opera  t  :  on .-71  Type  Ft  "nd  1  es  ,  WSLC  2 e  s  o  a rcF  V.bn.or iiiuuii: 

Kb .  lir;  Weapon  Systems  Evaluation  Group,  The  Pentagon, 
Washington  DC,  15  October  1939.  (AD  261  732) 

5.  Bridgman,  C.J.  Class  lecture,  27  May  IPSO.  Wright -■ 
Patterson  APB  Oil:  School  of  Engineering,  Air  Force 
Institute  of  Technology. 

6.  Droll,  S.D.  and  F.  vonllippcl.  "Limited  Nuclear  War," 
Scientific  American,  233:  27-36  (November  1976). 

7.  SACM  105-2,  Vo  1.  II.  Cl i mat ol og i c a  1  Wind  Factors. 

Omaha  NH:  Strategic  Air  Command ,  January  i960. 

8.  Glasstone,  Samuel  and  Philip  J.  Dolan.  The  Effects  of 
Nuclear  Weapons,  U.S.  Government  Printing  Office , 

W  a  sTfi  n  g  t  oiTb  •  L  ,  1 1 ,  i  r  d  Fdition,  1  97  7  . 

9.  - Handbook  of  Mathematical  Functions,  edited  by 

Milton  Abramov: iT t  and  "l  reno*  Jtcgun  I: .  S.  Government 
Printing  Office,  Washington  DC,  June  1964. 


61 


Appendix  A 


Fort  ra  i:  do  V  A' '  !.T  ! 

Tlij.?  appendix  contains  t  he  I  r >  1  lowing  items: 

--  a  copy  of  the  code  MULTI  which  computes 
fallout  contours,  and  a  glossary  of  the 
terms  used  in  MULTI; 

-  -  A  user's  guide  to  running  the  compiled  version 
of  the  MULTI  program; 

--  a  sample  of  output  from  MIU.TI  and  an  explan¬ 
ation  of  how  to  interpret  the  output. 


MULTI  Code 

MULTI  is  a  derivative  of  SMLA'd,  a  single- burst  fallout: 
code  written  by  Dr.  Charles  J.  Bridgman  of  the  Air  force 
Institute  of  Technology.  Several  changes  were  incorporated 
into  SMEAR  to  produce  MULTI,  most  notably  the  multiburst 
f(y)  distribution  developed  in  Chapter  II.  MULTI  was  also 
designed  to  be  run  interactively,  with  inputs  coming  from 
both  the  user  and  an  attached  data  tape  (TAPF.4).  The  user 
inputs  are  the  parameters  of  the  particular  problem  to  be 
studied,  while  the  data  tape  consists  of  a  table  of  coef¬ 
ficients  supplied  by  Colarco  (Ref  5:67-71).  The  table  of 
coefficients  is  required  in  the  AFIT  calculation  of  g(t)  . 

The  operation  of  MULTI  is  st  ra  i  glit  forward  ,  with  all 
dose/ do so  rate  quantities  being  calculated  as  described  in 
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Chapter  11.  The  contour  coordinates  are  given  in  terms  of 
x  (downwind  distance  from  line  of  hursts)  and  y  (cross- 
wind  distance  from  center  of  field).  These  contour  coordi¬ 
nates  are  displayed  with  t  lie  output  at  the  terminal,  and  are 
also  written  on  a  data  tape  (TAPF.6).  This  data  tape  can 
then  be  used  as  the  input  for  contour  drawing  routines. 

The  generation  of  contour  coordinates  is  an  iterative 
procedure.  The  downwind  distance  x  is  increased  by  a  small 
amount,  and  the  dose  (or  dose  rate)  at  that  x  coordinate 
and  the  center  of  the  field  (or  y  -  0)  is  then  computed 
and  compared  to  the  desired  dose.  If  t lie  computed  dose  is 
equal  to  or  less  than  this  desired  dose,  we  are  either  on 
or  outside  cf  the  desired  dose  contour  and  the  y  coordinate 
is  assigned  a  value  of  zero.  If  the  computed  dose  is  greater 
than  the  desired  dose,  we  are  within  the  desired  contour  and 
the  y  coordinate  is  increased  by  a  small  amount.  A  new 
dose  is  computed  for  this  new  (x,y)  position.  If  this 
computed  dose  is  still  greater  than  the  desired  dose,  the 
y  coordinate  is  again  increased.  This  procedure  continues 
until  the  computed  dose  equals  the  desired  dose;  1  lie  y 
coordinate  is  then  the  total  distance  incremental!;,  t-avclcd 
from  the  center  of  the  field  to  this  equality  point.  The 
downwind  distance  x  is  increased  again,  and  a  new  y 
coordinate  is  determined  using  the  same  procedure.  Twenty- 
six  pairs  of  (x,y)  coordinates  are  generated  to  define  a 
desired  dose  contour. 
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A  glossary  of  terms  used  within  the  code  is  provided 
in  Table  A-l.  These  terms  are  presented  in  order  of 
appearance  within  the  program. 

User's  Guide  to  MULTI 

Table  A-2  shows  how  the  parameters  of  a  problem  arc 
entered.  The  first  directive  sets  up  the  internal  switches 
for  the  program: 

M  determines  the  type  of  g(t)  calculation 
(0  =  1\’SI:G ;  1  =  A  FIT)  ; 

--  MD  determines  the  type  of  output  desired 
(0  =  dose  rate;  1  =  dose); 

N  determines  the  number  of  contours  to  be 
generated. 

The  parameters  of  the  scenario  are  entered  next,  with 
the  desired  contour  values  entered  first.  The  yield  of  the 
attaching  weapons  (in  megatons)  and  their  fission  fraction 
are  then  requested,  as  are  the  average  v/ind  velocity  (in 
miles  per  hour)  and  crosswind  shear.  The  number  of  bursts 
must  be  entered  next,  along  with  the  crosswind  width  of 
the  field  (in  miles). 

If  the  desired  contours  are  to  be  computed  in  terms  uf 
dose  (that  is,  MD  =  1),  the  total  number  of  hours  'Her  which 
the  dose  is  to  be  integrated  is  entered  next.  ih.v  integra¬ 
tion  over  time  is  the  Way-Wigner  approximation  ,!Vi  S :  .194) . 

If  a  time  of  infinity  is  desired,  zero  is  the  regained  input. 
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Finally,  if  all  calculations  are  to  be  made  using  the 
AFIT  g  (t )  (that  is,  M  =  1),  the  parameters  of  the  activity- 
size  distribution  curve  must  be  entered.  To  approximate  the 
DF.I.FIC  default,  the  parameters  given  in  Chapter  111  should 
be  used. 

Interpretation  of  Output  from  MULTI 

Table  A-3  is  the  output  of  the  scenario  entered  in 
Table  A- 2.  A  short  summary  of  some  of  the  input  data  is 
first  given.  "AFIT  CALCULATION"  means  the  program  will 
utilize  the  AFIT  g(t)  in  its  computations.  The  yield, 
fission  fraction,  wind  and  wind  shear  are  printed  out,  as 
is  the  desired  dose  contour.  The  26  pairs  of  coordinates 
for  the  contour  are  then  listed.  Note  that  this  listing 
is  just  "half"  a  contour;  due  to  the  symmetry  of  the  cloud, 
the  other  "half"  of  the  contour  has  the  same  x  coordinates 
but  the  negative  of  the  y  coordinates. 


ooooo  o  oo  o  ooooon  o  o 


c:oMi’im:K  com: 


pe  :> o r s  v  vul t i  ( ^  - a,  t^p.4) 

"  m.  IS  ■  *  T  :  '  .  -  .  ■  •  •  :  D  r-  71  OVD. 

TAPES  IS  At-  OUT  3 1'  1  TA-S  rDR  CONFUTED  VALJES. 

oi  4-enei  r\'  x  c 2 r ) ,  Ym 

C0uMCM  D(«),jC(;  )  ,  D  (  '  > ,  A  L  ,  n'r 

CCMMCN  T  3,  SI  3  0, 5  I  CH,  SH'  ,S‘C,  V  X ,  7  n,7,“  ,TP  ,K,  M,MR 
COMMON  NN,w,TD 

INTERNAL  SWITCHES  = 

H=:  FOR  'J>  EG*  >•-  1  FOR  a  FIT 
M0=t:  F  OF  33SE  k  i  T  =  C»/H'-  ),  k’0-l  rpR  003  E(R> 

M  IS  THE  NU  ' 3 E  R  OF  COiTDLF  <  TO  BE  COMPUTED 

PR  I  N  T  *  fMi<!TER  Hi  «-  OR  1)  ,MDU  OR  1 )  »  A  NO  N<  «  CDNTOJRSm 
RE  AO  '  ,  M,  MO,  M 

RE  A  0  N  VALUES  OF  CC  MODES  ;rOUEST  ED  IN  R/HR  OR  R 
PRINT*  , ••ENTER  CONTOUR  7ALUES=  ** 

REA  0* ,  (C(I) ,1=1, N) 

READ  Y  2  F.L  n  (IN  *7), FRACTION  OF  YIELD 
WIND  (IN  HI'-IR)  ,  AND  WIND  SHE  A  (IN 
PRINT'  j  **  E  'T  T  E  f  y:  elf,  fra  CTION,  WiNC 
RE t  0  *  f  YM,  FF  ,  VX  ,SWR 
PRINT*  ,**EU  ER  i>  CF  DIJRSTS=  " 

READ'-  ,  NN 

PRINT1  >  **  H  T  £R  N-S  WIDTH  OF  FIlLC  = 

RFC  O'  ,  K 
W-W/2. 

IF(MP.  £0.1  . )  pfihT'V'ENTE.p  TIME  OF 
IF<HD.  EC.i.)Rf/-?i,T-' 

SC  (SOURCE  NORMA  LI T AT: ON  CONSTANT)  YIELDS  R  OR  R/HR 
IF  (  M*  EO  ,  •)  S  C  =  2 »  <E  ♦  ’  V-fp-'v^ 

IF<M.t  D.i.  )  SO  =  2.  '5* £+••-'  *FF*  Y M 

YIELD  DEPEND  mt  DON:  'ART  S  '  r-  NOW  CO'=U''ED 
HC  IE  IN  < I L  0 f E  T 
SIGH  m,r  SI  GO  t,=  I  N  GTU'JTE 
TO  IS  IN  HOURS 

HC=  *•  LLnG(Y»-)  -.  *»  '  rUl  Or-(YH\ +B  .  S  )  ■*  -.  n  S  (  A  _  'j  G  (  /  R  / 

<  *  2 « -  2 ) 

SIGH  =  .  1  !'  HC/f  .  ?  :: 

SIr-  0  =  -Vm.r  +A  .  C  Drv-M  / -»,  -  •  .  S  l  /(■••-*■  (A  L  C  G  f  Y 11  >  *■  > .  r  '  '  *  )  ' 

C  TC  =  1  .  '  V?*  (12.  -T/  . L  PC/*  i  -  (1  ‘  X-  t 

K  -PC ' H'/E2:  .)  ) 

C  SET  A.PHA  A  NO  BETA  F(p  w5r.  G  CALCULATIONS 
AL  -  ■»  A  . 

BT  =  .09 

C  IF  THIS  IS  G  HS'G  C  ALC'IH-  ICC,  S<I 3  TO  £ . 

IF  ( M, EC. ? )  SOTS  6 

PRINT*,  "FMEP  Alp  PA  AND  BETA  =  •• 

RE  AO  *,  AL,BT 


AS  FISSION, 

HR-i  ) 

level,  and  wind  i H T  = 


OUKATIONs  ** 
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o  o  o  o  o 


C  OOWCFT  <  I  _  0  F  -  "T  TD  *  7 L  0 v  E  T  c  R  S 

7  V  -  4  ~  '  a/t  - 

C  DERIVE  DONSTAMS  FQ?  C(T) 

DO  ?S  J-2,b 
C  (  J-  1 )  -  COT  FF(7K,  J  ) 

25  CONTINUE 

DO  25  J=-,9 
C  (  )-2)  =CC5  FFPK,  J  ) 

26  CONTINUE 
PRIM  lb 


■c 


c 

c 

c 

c 

c 

c 

c 

c 


55  FORMAT  H  ;  f:t  z\  LCULATTCN  ♦) 

Ip  TMIS  IS  A  '•SEC-  C  ALOLAT  ION,  PRINT  THAT 
ec  if  <m« en.: i  printed 
65  FORMAT  (  *  4 ESC  CALCULATION  •) 


PRINT  CUT  TH-  IN-f-U*  PS~Ali:cRS 
PR!  N T  ?  .  ,Y4,  f- 


70 

FQ:  -  AT 

( 

YIELD  -  ,".2, 

M r  . 

FIS'  I  O'!  r 

'•.AO  T  I  ON-  =  •,-3.2) 

PRIM 

7  •: 

> 

vx, S4r 

75 

F0.-MA1 

(  ♦ 

N[  NO  =  *  ,F-.  i,  ■ 

WIND  5  He-. 

=  4  tF>  .?»  ’ 

FP.I  •>.  T 

t  i 

8C 

FO-  ‘Ml 

(* 

\  D  r I  VI  i  Y-~  T  7  E 

r  r;  ; 

I  r  U  T  i  0  '  i  ~ 

lO  GNO'IH  L  '1 

POT  = 

1  * 

9 

PR*  'J T 

2 » 

T 

82 

FORMAT 

( i  y 

v.  A  :3  * 

N  A  _  r 

u;  =  »  ,  r  r  « 

1 AND  -  ~ET  A  *  * 

^  F-.J 

) 

FI  HD  THE  TI'-'S  0-  MAXIMUM  G  (T)  AMP  CALL  17  "TM” 

86  CALL  GMAy 

DO  12'-  K  =  i,N 

FIND  U°3L3=E  i  MESSED  1  ION  EY  CALLING  "UPHILL" 
CALL  JPHILL 

RETURN  •*  I  T  H  VtLU-  FD  =  T  D 
IF  (13.50.--  1.)  G0T0  1 2  v 
FIND  no-(NSLCF-  INTERSECTION  3Y  VALUING  "DWHILL” 
CALL  DWHILL 

RETURN  -» I T H  vc  LU-  FDR  TF 
DIVIDE  HOT  LI ‘41  INTD  2*  INURVALS 
OT  =(TF-TE) ‘ .  4 


) 


PRINT  T  A 3  L  3  TITLE  FOr  APPROPRIATE  CALCULATION 
IF  HD.  EC.  )  p^INT  ?:  ,  0(0 

90  FOF  H  A  T  (*  COORDINATES  FO  fc,F7.1,*  R/  HR  IS  ORATE  .  I  N  E  *  ) 

IF  ( HD . N  E . .  )  °  R I N '  35,  0 (<) 

95  FOPMAM-  C  0  ORCI NA  T  rS  F0‘  ',r7.1,*  &  13^^'  LIVE4) 

DOPPUT-:  AND  p;i*r  X,x  VAlUv.1  FOR  THIS  CCiM'UR 
X  IS  ODWNWirO  01 S P A N Cr  FRON  5URSTS 
Y  IS  C  ?  D  S  S  W 1  v D  DISTANCE  FROM  UDT  LINE 
DC  131  I  =  1,2'. 

T  =  FL  CAT  (t_l  )*oT  ♦  T3 
X(I)  =  VY*r 

C  rlND  COSE  ON  THE  HOT  LINE  rCR  EA  C 4  ST  E=>  X  (I) 

Y  ( I )  = 

7  =  A  BE  (  CY(T)-M)  /  SIGY  (X(I)  )  ) 
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CUM.  NORMAL  FUNCTION  AS  GIVEN  nY 
ACRti-'o^ir7  anO'tegjm  ?.;?2 

F  =  1  .  /  (  2.M1.  +  .19!  5d  V7J.  .  1  ibiD^  74  I 

*  f.C  .•  i'3k'.*7**?+,.i9927*7*,:-)-i'*A) 

ra=i.-F 

rF=FA-F 

FY=NN/  (2k  ,-0  >(r  ° 

D=SC'G(T)/'K*  F Y 
I F  (  HP,  NE.  ’)  p  =  C*E 1 <  T) 

IF  QOS  "  3  N  THE  *3  7  LINE  ]S  GREATER  THAN  THE  DESIRED  DOS- 
GO  CROSSWINP  UNTIL  DESIRED  DOSE  IS  I NT  ERSEO  TED. 

IF  (P.  GT.O  (<)  )  GOTO  U, 

/(I)  =  P. 

GOTO  15 

D  ETERMI NA  T I  ON  C-F  CRDSGWI'D  DISTANCE  NEEDED 
TO  Ill  TEE  SECT  OESI  RE  n  DOSE 

id"  a=g 

MA)(W=(  W+3*  STFVfXd)  )  )  *  13  . 

DO  It;  J=t,MAXW 

6=A+.i 

?  =  ABS.  (  (A-M)  /  STGY(X'd))) 

F  =  1./(2,B(1,  +  .i9t£>-*'7+.lld5^2*2 

*  +.  31  •  3  VV  7'*  7  T.C  !  SI  c'l  *  7*- .'»)*-»  4) 

IF(A.LT.W) F=  t • -F 

FY  =  NFV  <2"  M)  »- 
0  =  50" G  (T) /7X* 

IF<  MO,  NE,  r  .>  0=0*  BUT) 

IF(0«GT  .  j  ( K)  )  GO  TO  lv6 
Y  ( I )  =  A 
GO  TO  i T- 


106 

CONTINUE 

iC5 

PRINT  lt',I,X  (1)  f  Y  ( 1) 

WRITE (  f » 1 1 2)  ,  X  (I)  ,  Y  (I) 

1  i'¬ 

FORMAT  (I>f2Pa.l) 

ll? 

FORMAT  Ci'X,2-  j‘,2) 

Ill 

CONTINUE 

PRINT  HE  FIML  POINT 
X  ( ? G >  =  VX*  Tr 
Y(2M  =  C  . 

PRIM  1U  ,?F,Y(23)  ,Y(EO) 
M^ITE (1 ,  11?  J  ,X<  25)  ,v  <26) 
123  CONTINUE 

STOP 
END 


J 


SUB  RCJT  1  ME  3MAX 

cc-'^c-;  r  (♦),:«(-  d  ,  or  )  ,al  ,'  * 

c.C:  !-l\  tc,  hi  :-o, :  :gu, sm;  ,sc,  •/x,t:.,tu,tf,<,  m,.m? 

COMMON  K\’»  W ,  T  D 

C  MAXIMUM  G  <  T  )  15  FOUND  QY  T  FACING  GC“) 

C  TO  IIS  PEH  I  SI  .1  HOUR  INCREMENTS 
DC  ♦  J=i,i  I 
T  =  r _  3 AT  ( J- 1  )*  •  1 
OH(  J)  =3 (Tl 
IF  ( J.E  0.1 )  GOTO  A 
IF  (DM ( J) . LT.OM ( J -1) )  GCTG  5 

4  CONTINUE 

5  TM=T-.i 
RETURN 
END 

Q  +  +  +  *¥**•■  •  **  +  +  *■:  *UI4»4H'  *++++*  +  *  +  +*+■(■’'■* 

SUBRCJTI  !E  J°HILL 

COMMON  0(4),  CMC  '))  ,  C  (/  )  , A  L , FT 

COMMON  TC,SrSG,  SIGH,  SHF  ,SC ,  VX  ,  TE ,  TM  ,  TF ,  K,  M ,  MO 

COMMON  NN, W, T  0 

C  THIS  5  UrROUT  THE  FINOS  THE  CQIM  CN  THE  HOTLINE 
C  CLOSEST  TC  THE  BURSTS  WH'-^S  THE  DESIRED  CONTOUR 
C  OOSE  INTERSECTS. 

C 

C  TB  =  -9S.  I?  A  “LAG  TO  INDICATE  THAT  J=1D  IS  LESS  THAN  T  M 
TB  =  -9?:. 

C  T0=INITIAL  TIME 

TO  =  -(3,»SIG0)/VX 
IF  (M.NE.T)  TO  =  0. 

T=TQ 

C  SET  OT  FOR  N3EG  (FDA  T  <  7 ‘SIGMA ) 

DT  =  .1*  A  3  5 ( T  0 ) 

0*(1>  =  . 

DO  2  J  =  2 , 1.  •’ 

C  SET  OT  r3R  A  FI ”  0°  RESET  IT  FOR  WSEG 

IFCf.3E.AP5  (TO)  )  0?=,  1 
T  =  r  +  OT 
x  =  r»vx 
?=W/SI3Y(X> 

F  =  1.M2.  ‘  C  ,*.1053  c'.*  Z+.  11E  IDG  “Z*  Z 
*  +  .C;)344«L?'1*3+.Jl?r??''7'*4)  «1G) 

FA=1. -F 
FB=FA-F 

FY=SN/  ( 2  •  '  W  ) *f ^ 

OM  (  J)  =  S3  *  G  (  T )  /  V X *  F  Y 
IFCMI.NE.II  Du(  !)  =  PI(T)*CM(J) 

IF  OM(J)  .GE.D(K)  )  GOT  0  3 
IF (7.GT.TM)  RFTURM 

2  CONTINUE 
R“TURN 

C  INTERPOLATE 

3  DLT  =  DT*CDM( J)-D(X) )/(  CM ( J) -DM ( J-i  ) ) 

TB=T -DL  T 

RETURN 

END 


oooooo  o  o  ooo 


£  *  •  *  #**♦.•<  1*»  i**  -  M  *>  •»  )  '  »ue«s  4AA.»*. 


Sinc-cuTI  .r  jJHTL1. 

CO“MCN  LM»),jh(;  0  , C<")  ,  A.  ,  ■  T 

COMMON  TO,  51  30,  S  IGW,  SHR  ,~C,  7X,Tc,TM  ,TF,K,  M,  MD 

com mon 

THIS  SUBROUTINE  PINTS  THE  c?IH  CN  THE;  HOTLINE! 

FARTHEST  FROM  THf  njE-TS  WHERE  THE  DESIRED  DOGE  I  N 1  E  R  S 


JL  =  TM/.5  ♦  2. 

DO  *  J  =  JL,^ 

T  =  “LOAT  <  J-l)  *.E 
x  =  r*vx 

7=HXSI GY (X ) 

F=l.  /  (2  .M  l.+  .  19f ,F,EA*  7+  ,  he  19  '<*7+2 
*  +  .  )  3  3  4<<  *  7  -  *2  +  .  •*  i  9  -  2  7  *  7  '  <  k  )  *  ■»  4  ) 

FA=l.-F 
FB  =  F  A -F 
FY=4N/ ( 2. * W)« 

DM ( J) =SC*  G  (T) /VX‘  FY 
IFMD.SE.:)  0^  (  J)  =nI  ( T) *  Dw ( J) 

IFO(K)  .GT  (  J)  >  GOTO  G 

4  continue 

INTERPOLATE 

5  DLT  =  .F- <0 (<> -0 M( J)  )/ (  Cu( J-i ) -DM ( J) ) 

TF  =  T  -DLT 

RETURN 

END 

FUNCTION  3(5) 

COMMON  PC*),  DM  (5  0)  ,  0 ( 7  ) ,  AL  ,  ?T 

COMMON  TO,  5130, SIGH,  3H-  ,SC,  VX  ,TE  ,TM  ,TF  ,.K,  M  ,  MD 

COMMON  NN,W,TD 

THIS  SUBROUTINE  CC-FUTFS  G  <T)  FOF  ATIT  OR  WS  EG . 

IF  (K.EO.:)  GOTO  R 
IF  (S.LT.,1)  GCTC  3 

AF  IT  DOKFUTUICNS 

R  IS  PARTICLE  RAOIJS  (TN  METE-S) 

R=C  (i)  /  (3»  *  S)  *C(  2)  /(  3*  ~  4)  *0  (3)  /<S**  3) 

R=9+C(R> /(?*  S) +:  (5) /SfC  (C ) +  C<? ) /SORT (S ) 


C  CONVERT  METERS  TO  MICRONS 
K  =  RM.  i 
C 

C  A  IS  ACTIVITY-SI 7E  FJMCTIDK 
A  =  t  . 

P  =  ( A L  C G ( ; )  -  A 10G(AL)  )/9T 

A  =  EXP(-.  5*P»F)  /(SORT  (  G. 2;  3) J  GT*R) 
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c  or/dt  is  DEsmrivs  of  r 

D?OT:-:.  ,  +  Z  (1)  /  (3  V. • )  -  4.  *C  (2)  /  (5*  r 5)  -  3.  *C(  3)  /IS 
0R0T  =  0RrT-2  .  r-C(0/(S*-‘')-C(:')/(S*S)-.:^C(7)/(;“l.?> 
DPDT  =  C'  OT*  1 .  F  *■ <  6 
G=  A'ABSORDT) 

RETURN 

C 

3  G— .  • 

RETURN 

WS  EG  Z  OHPUT  ATI  ON’S 

4  X-VX‘S 
0=A6S(X) 

H=i  ,  +  .19-:.35k  '  Q+.  Ii3i9'4'  0*Q+«:  3  +  .  wi952  7*Q* F  t 

PHI  =  i,  -  l./(2.»'  h*K,0 
IF  (X.LT.j.)  PHI  =  1 • - c  H I 
G  =  PHI  *  SK=>  (-S/T  C)/TC 
RETURN 
EN’O 

*>■*».#■**  *►*****«.- 

FUNCTION  31 (S) 

COMMON  F(  4  >  ,OM(i  <.*}>  ,0(7  ),AL,PT 
COMMON  TO,  513  0, 315'',  SHE  , SC,  VX,T?,T>!  ,7F,K,  M,M3 
COMMON  NN,W,TO 

THIS  FUNCTION  CONVERTS  HOSE  RATE  TO  003E 
USING  THE  vny-WIGN-R  A  3SC  OXI  MAT  I  ON 

TP  =  S 

IF(TA.LT, . 1) T A  =  ,  1 
TE  =  T/WTO 
BI»F./(Ta**.2> 

IF(TE.GT.TA)  31  =  9  J-E./(T  Ev-*.  2) 

RETURN 
END 

Q*4«44444r4H<*  •  f  *  *•  *  ■»  ■*  »  ♦»**♦.***»»*««*».  »**#.4*v»*  ♦*#¥*** 

FUNCTION  3  I  GY ( X ) 

COMMON  0(4),0»'Ci  0  )  ,  ZC  ),  AL  ,?T 

COMMON  TC,S1G0, SIGH, S-l?  , SC,  V*,TC-  ,  rr  ,  .  M ,  MO 

COMMON  N  J »  d  ,  T  D 

C  THIS  rUNCT ION  OOMPJTES  SIS**  Y,  AS  F£E  WjtlO-l. 

C 

TS=X/VX 

IFCTS.5T  .3. >  TS*3  . 

TP=i .♦6.,rS/TC 

SIGY  =  Sn3T ( S I  GO '  *2*  TR  ♦  (SIGH*  SHR-*  */ V X ) *  *  2 ) 

RET  UR N 
END 


71 


o  o  o 


»*  *»>  M  M  *  t  f  ¥  t  »4».f  **»*.*♦  |  +  »<4M.*V 

FUNCTION  :otrF(H,N) 

DIMENSION  A  (2!r  :  ,  9) 

HIS  rUNCTIjv)  REACi-  TJC  DATA  TAPE  AND  INT ERPOL ATES 
FOR  THE  CCE-F1DIE  NTS  NEEDED  FOE  THE  GIVEN  ALTITUDE  *< 

REWIND  ft 

00  II  r  =  i  ,230 

RE  P  D(- ,1  ;  >  (ft  u,  J)  ,  J=i,  «•> 

RE  A  0  ( V  » 1  .  I )  (ft  (I,  J),  J  =  o,  o) 

1C -1  FORMAT  (F-  .1,4  Ell  .5) 
lQi  FO E  M ft  T ( F  . .  1 , 3 E 1 1  . F) 

IF(?<.LE.ft  (I,  i))  GO  TO  11 

10  CONTINUE 

11  X2=A(I,1) 

Xi=A(I-i,l> 

Y?  =  A(I,N> 

Y1  =  A  (1-1,0 

COEFF=( (Y2-Y1) /( X2-X1)  )  *7K  + ( ( X2* Y1 - Xi* Y 2)  / ( X2-X 1) J 

RETURN 

END 
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TAR Lh-  A-  1 


Gloss ary  of  Terms  V.ilhin  MU!T  1 
INPUT  PAR  AM  LTF.RS 


DCI) 

Contour  value  in  dose 

or  dose  rate 

YM 

Yield  of  the  weapons 

(in  megatons) 

FF 

Fraction  of  yield  due 

to  fission 

VX 

Velocity  of  the  wind 

(in  miles  per  hour) 

SHR 

Crosswind  shear  (in  hours  ^) 

NN 

Number  of  bursts 

W 

Width  of  the  field  (in  miles) 

TD 

Time  of  duration  for 

dose  computations 

(in  hours) 

AL 

Alpha  parameter  from 

act iv ity -  sice 

distribution  curve 

BT 

Beta  parameter  from  act ivity- s ice 

distribution  curve 

WSEG-10  PARAMETERS 


HC 

Nuclear  cloud  height 

(in  ki lof eet ) 

SIGH 

Cloud  thickness  parameter 

SIGO 

Initial  cloud  radius 

pa  ramet er 

TC 

Time  constant 

SC 

Source  normalization 

const  ant 
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TADI-i:  A  - 1  (font  ‘  A: ) 


roMi’UTAT i on  paranftnng 


C(J) 


MAXW 


Nuclear  cloud  height  (in  kilometers) 
Coefficients  for  AF1T  g (t )  computation 
Interval  along  hot  line 
Time  of  maximum  g (t ) 

First  time  of  occurrence  of  desired 


dose  on  hot  line 


Final  time  of  occurrence  of  desired 


dose  on  hot  line 


Time  (in  hours) 

Downwind  distance  (in  miles) 
Crosswind  distance  (in  miles) 
Cumulative  normal  function 
Crosswind  distribution  of  activity 


Dose  or  dose  rate 


Crosswind  mileage  counter 
Crosswind  limit  of  commit  at  ion 


FUNCTION  G  (S)  FAN.'.' ,TTFRS 


DRDT 


L 


Particle  radius  (in  meters) 
Activity-size  function 
Change  in  particle  radius  with  respect 
to  time 

Cumulative  normal  function 


TABU:  A- 2 


Interactive  Inputs 

ENTER  M  (0  or  1),  MD  (0  or  1),  AND  N  (#  CONTOURS)  =  1,1,1 
ENTER  CONTOUR  VALUES  =  15 00 

ENTER  YIELD,  FRACTION,  KIND  LEVEL,  AND  LIN'D  S!!R  =  1  ,.5,35. 
ENTER  #  OF  BURSTS  =  4600 

ENTER  N-S  WIDTH  OF  FIELD  =  ISO 

ENTER  TIME  OF  DURATION  =  0 
ENTER  ALPHA  AND  BETA  =  37,1.528 


75 


TAP.i  r  A- 3 


Sample  Outp  1 1 1 


Vita 


John  P.  Crandley,  Jr.  was  born  on  25  November  1951  in 
Pittsburgh,  Pennsylvania .  He  graduated  from  Milford  High 
School  in  Milford,  Connecticut  in  1969,  and  attended  the 
U.S.  Air  Force  Academy  where  he  received  the  degree  of 
Bachelor  of  Science  (Chemistry)  in  June  1973.  After  grad¬ 
uation,  lie  attended  Undergraduate  Navigator  Training  at 
Mather  Air  Force  Base,  California,  which  he  completed  in 
May  197-1.  He  then  attended  B-S2  Bombing/Navi  gat  ion  School. 
He  was  assigned  as  a  B-52G  Navigator  and  Radar  Navigator 
between  January  1975  and  August  1979  at  Griff iss  Air  Force 
Base,  New  York.  Fie  entered  the  Strategic  and  Tactical 
Sciences  program.  School  of  Engineering,  Air  Force  Institute 
of  Technology,  in  August  1978. 
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